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We develop an accurate finite-time scaling analysis of the angular width of the coherent backscat¬ 
tering (CBS) peak for waves propagating in 3D random media. Applying this method to ultracold 
atoms in optical speckle potentials, we show how to determine both the mobility edge and the criti¬ 
cal exponent of the Anderson transition from the temporal behavior of the CBS width. Our method 
could be used in experiments to fully characterize the 3D Anderson transition. 
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In disordered media, the absence of diffusion arising 
from the spatial localization of single-particle states is 
known as Anderson localization (AL) [I,]. In three dimen¬ 
sions, AL manifests itself as a phase transition, which oc¬ 
curs at a critical energy, the mobility edge (ME), separat¬ 
ing a metallic phase where states are spatially extended, 
from an insulating one where states are localized. The¬ 
oretically, much efforts have been devoted to the study 
of the critical properties of the Anderson transition, such 
as wave functions at the ME m or critical exponents 
[4| . In practice however, only a handful of experiments 
have found evidence for the three-dimensional (3D) An¬ 
derson transition [M3- For matter waves, its critical 
(universal) features have been only investigated in the 
context of quantum-chaotic dynamical localization Q 
but no such experiment in 3D disordered potentials has 
been reported to date. 

In addition to the intrinsic difficulty of achieving wave 
localization in three dimensions, one reason for the 
rareness of experimental characterizations of the Ander¬ 
son transition lies in the lack of easily measurable ob¬ 
servables displaying criticality. In the context of atom 
optics, a routinely used approach consists in tracing the 
evolution in time of the spatial width of a spreading 
wave packet 0-E3- While AL implies a saturation of 
the width, the contrary is not true as classical effects 
can as well entail a saturation or a slowing down in time 
|T3 . Furthermore, atomic wave packets have rather large 
energy distributions even when cooled down to very low 
temperatures, which forbids an accurate resolution of the 
critical region around the ME. Thus, any exploration of 
the Anderson transition with cold atoms should ideally 
be complemented with a clear demonstration of phase 
coherence, and should achieve a good energy resolution. 
For the latter issue, a first step has been reached in recent 
measurements of the ME based on a frequency modula¬ 
tion of the disorder Q although the experimentally 
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measured ME seems significantly higher than the one 
predicted from extensive numerical calculations E1E1- 
and upcoming experiments are moving toward a genuine 
filtering of the energy distribution, required to access the 
critical properties of the transition [15[ • For additionally 
proving genuine phase-coherent scattering, the coherent 
backscattering (CBS) effect is a promising tool. CBS has 
already been observed in several experiments with cold 
atoms 0 , light [ 13 , Ell, acoustic 0 or seismic waves 
[23 in the (metallic) regime of diffusive transport. Inter¬ 
estingly however, CBS shows up not only in the metallic 
phase, but all the way across the Anderson transition. 
The question then naturally arises whether the CBS peak 
itself could be used as an observable for accessing the crit¬ 
ical properties of this transition, in which case one would 
simultaneously ensure phase coherence. 



FIG. 1. (color online) Contour plot of the averaged momen¬ 
tum distribution of a matter wave, obtained after propaga¬ 
tion of a plane wave |fco) (fco = 0.6ea,) in a speckle poten¬ 
tial of strength Vo = 1 for a duration t = 800. The propa¬ 
gated state is here filtered around energy E = —0.4 (metallic 
regime). The CBS peak, of angular width 2A#cbs, is visible 
at fc = — fco. Here momenta, energies and times are respec¬ 
tively in units of £ -1 , h 2 /(m(,' 2 ) and m( 2 /h, where £ is the 
correlation length of the potential. 


In continuous-wave optical experiments, it is known 
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that the CBS lineshape changes at the critical point [2lj . 
Unfortunately, this feature is usually smoothed by ab¬ 
sorption or finite-size effects and cannot be used in prac¬ 
tice. In this letter, we explore the dynamics of the CBS 
effect in momentum space - in contrast with usual setups 
that search in configuration space - around the Ander¬ 
son transition. By scrutinizing the dynamics of the CBS 
angular width, A6 ?cbSj hr combination with a numerical 
filter that provides a high energy resolution, we demon¬ 
strate that A$cbs can be used to characterize the critical 
properties of the Anderson transition. By developing an 
accurate finite-time scaling analysis of the CBS data, we 
verify the one-parameter scaling theory of localization 
locate precisely the ME and extract the critical ex¬ 
ponent of the transition. We determine these parameters 
for a speckle potential, and find a good agreement with 
the predictions of the transfer-matrix method. 

As shown in j23|, CBS of cold atoms can be observed 
by tracing the evolution a quasi-plane matter wave in 
momentum space, a proposal recently realized experi¬ 
mentally [16;]. Let us thus consider a matter wave ini¬ 
tially prepared in the plane-wave state \ip(t = 0)} = |fco), 
and subjected to a 3D random potential V(r). Follow¬ 
ing experiments, we choose V(r) to have the statistical 
properties of a blue-detuned speckle pattern. It is cus¬ 
tomary to shift all energies by the average value Vo > 0 
of the speckle potential, leading to the on-site distri¬ 
bution P(V ) = exp[-(U + Vb)/V r 0 ]B(U + Vo)/V 0 (0 is 
the Heaviside function), and the two-point correlation 
function V(r)V(r') = U 0 2 [sin(|r — r’\/(/)/(\r — r'\/C,)] 2 , 
where £ is the correlation length. In order to accurately 
pinpoint the ME E c , it is useful to restrict the evolu¬ 
tion to a narrow energy range (±er) centered at a given 
value E that we wish to tune around the E c , by apply¬ 
ing a Gaussian filter [24| exp[— (H — E) 2 /(2a 2 )] (where 
H = p 2 /(2m) + U(r)) on the initial state |fco)- This filter 
makes it possible to accurately extract E c which other¬ 
wise would be smoothed by the natural energy distribu¬ 
tion of the initial plane wave in presence of the disordered 
potential |25|. Throughout this Letter, lengths, mo¬ 
menta, energies and times are given in units of £, £ -1 > 
h 2 /(m( 2 ) and m( 2 /h, respectively. We discretize the 
Hamiltonian H on a 3D grid of total volume (60 x tt() 3 
with periodic boundary conditions. Each cell of size 7 t£ 
is divided into 2 steps in all three directions. In the fol¬ 
lowing, we use Vo = 1, a = 0.02, fc 0 = 0.6. 

The temporal evolution and the filtering are performed 
using a Chebyshev scheme. The evolution operator over 
At, e ~ lHAt / h [resp. the filtering operator] can be ex¬ 
panded in a series of Chebyshev polynomials of the first 
kind of aH+b [resp. a(H—E) 2 +b] with a , b conveniently 
chosen parameters - see [2a [27| for details - whose coef¬ 
ficients are Bessel [resp. modified Bessel] functions of 
argument proportional to At. The temporal evolution 
can be computed by iterating small time steps, each 
involving a limited number of terms in the Chebyshev 
expansion. The momentum wavefunction is obtained 
by Fourier transforming the final wave function \ip(t)). 


The procedure is repeated over 6 x 10 3 configurations 
of V(r ), yielding the averaged momentum distribution 
n(k,t ) = |(fc | ii>(t)) | 2 . We show in Fig. [T]the numerical 
distribution n(k,t) obtained at long times for an energy 
E = —0.4 which lies in the metallic regime E > E c . 
n(k , t) clearly displays a narrow interference peak of an¬ 
gular width A$cbs an d centered at k = — k 0 (in red 
in Fig. [U. This CBS peak sits on the top of a time- 
independent isotropic background (in blue in Fig. [T]), 
which in three dimensions has the shape of a spherical 
shell as a result of elastic multiple scattering off the ran¬ 
dom potential [23]. 



FIG. 2. (color online) Dynamics of the CBS peak across the 
Anderson transition. Left panel: angular width A#cbs versus 
time, in the metallic regime E = —0.4 > E c (green points), 
at the mobility edge E = E c ~ —0.48 (red points), and in 
the insulating regime E = —0.56 < E c (blue points). Right 
panels: cut along k x of the normalized CBS profile at three 
different energies. For each energy, profiles at three different 
times t = 2000, 4000 and 8000 are displayed, shifted with 
respect to each other for clarity. The CBS width rapidly sat¬ 
urates in the insulating regime, while it shrinks in time in the 
metallic and critical regimes. We find an excellent agreement 
with the temporal dependences predicted by Eq. ©■ 


We now study the time dependence of the CBS an¬ 
gular width, A0 cbs- Qualitatively, CBS is an interfer¬ 
ence effect between two waves that propagate along an 
identical multiple scattering sequence r±... r n but in 
opposite directions (28|. The interference term between 
these paths is proportional to cos[(fco + k ) ■ ( r n — tt)]. 
Therefore, denoting by Ad the angle (assumed small) be¬ 
tween k and — fe 0 , we infer that an interference is visible 
on average provided k 0 A9Ar(t) -C 1, where Ar(t) = 
(| rN(t) — ril 2 ) 1 ' 2 . We thus estimate the angular width 
of the CBS at a given time t to be A$cbs ~ l/[ koAr(t )]. 
The average distance between the first and last points of 
the scattering sequence depends on the nature of trans¬ 
port in the system. In the metallic regime E > E c: 
A?’(t) oc \JD(E)t with D(E) the diffusion coefficient 
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at energy E , while A r(t) oc t 1 / 3 at E = E c [29j and 
A r(t) cx £(E), the localization length, in the insulating 
regime E < E c . We thus have: 

(l/y/D(E)i E>E C 

k 0 A0 C bs - < 1A 1/3 E = E c (1) 

[l/£(E) £<E C . 

The time dependence of A$cbs is thus qualitatively dif¬ 
ferent in the three regimes of transport. In particular, a 
sub-diffusive behavior of the CBS width marks the posi¬ 
tion of the ME E c . We have performed numerical simula¬ 
tions of the momentum distribution for various energies 
E around E c ~ —0.48. We show in the left panel of 
Fig. [2] the CBS width as a function of time, for three 
different energies around E c . At long times, the results 
follow very well the predictions of Eq. ©• For each 
energy, we have obtained A(9cbs by first removing the 
isotropic background [24| from the 3D momentum distri¬ 
bution, then fitting the resulting momentum profile with 
a/[l + (fc + fco ) 2 //?] 7 (where a, /? and 7 are time- and 
energy-dependent fit parameters), and finally taking the 
half width at half maximum of the fitting function. Error 
bars on A(?cbs have been estimated from the standard 
deviations of /3 and 7 . We show examples of CBS profiles 
and the corresponding fits in the right panel of Fig. [2] 



FIG. 3. (color online) Left panel: diffusion coefficient (in the 
metallic regime E > E c ), and inverse of the localization length 
1 /£(E) (in the insulating regime E < E c ) versus energy E. 
Red dots are obtained from an analysis of the CBS width, blue 
squares from the transfer-matrix method, and green diamonds 
from the spatial spreading of a wave packet. Right panel: 
scaling function A versus energy around the ME, for various 
times ranging from t = 1280 to t = 7680. The curves cross 
at a common point E c ~ —0.48, which signals the location 
of the ME. Points are the results of numerical simulations of 
CBS, while solid curves are hts of these data using Eq. ©• 


According to Eq. ©, the CBS width is also propor¬ 
tional to the square root of the diffusion coefficient D(E) 
in the metallic regime, and to the inverse of the local¬ 
ization length 1 /£(E) in the insulating regime, which 
suggests an original way of measuring these quantities 
experimentally. To demonstrate the efficiency of such an 
approach, we have extracted D{E ) and 1/£(E) from the 


numerical data for A(9 cbSj by extrapolating the quan¬ 
tities l/[(k 0 A9 C Bs) 2 t] (for E > E c ) and l/(fc 0 A 6 > C Bs) 
(for E < E c ) to infinite times. The results are shown 
as red dots in the left panel of Fig. [Sj for various ener¬ 
gies around E c (no values too close to E c are shown due 
to the lack of accuracy of the extrapolation procedure 
at these energies. The vicinity of E c deserves a special 
analysis that will be described below). We have com¬ 
puted these quantities using the transfer-matrix method 
(blue squares) [30]. In the metallic region, we have also 
computed D(E) by yet another method that consists in 
analyzing the spatial width of a spreading, initially nar¬ 
row wave packet as a function of time (green diamonds). 
All the results for D(E) are in very good agreement. The 
predictions below E c tend to deviate far from the ME, 
which we explain by the difference in the definition of 
£,(E) in the two methods: the localization length that 
appears in A$cbs controls the exponential decay of the 
average density, whereas the localization length that ap¬ 
pears in transfer matrices controls the exponential decay 
of the average of the logarithm of the transmission [30]. 

Let us now explore the behavior of AOqbs in the close 
vicinity of E c . In this region, D(E) oc \E — E c \ s and 
£(E) oc |E — £’ c | _! ', where the two critical exponents v 
and s turn out to be equal for the Anderson transition 
in dimension 3 |3l|. Near E c , the three scaling laws © 
can be recast under the unified form 


A = 


1 

Lk 0 A9cBS 


F 


Xr{E)LV-‘ 


( 2 ) 


where Xr(E) oc E — E c , L = [t/(2nhp(E))] 1 ^ 3 with p(E) 
the density of states per unit volume at energy E, and 
F is a function characteristic of the transition. Although 
the system a priori depends on two parameters E and i, 
Eq. © thus suggests that A is in fact a function of a 
single parameter, and is therefore a good candidate for 
developing a single-parameter scaling description of the 
Anderson transition [22j. The introduction of the length 
scale L [32f allows us to establish a straightforward anal¬ 
ogy with the usual scaling theory of Anderson localiza¬ 
tion for time-independent disordered systems d, [22j . A 
direct consequence of Eq. © is that when In A is plotted 
against E, the curves at different times should cross at 
E = E c . This behavior is well visible in the right panel 
of Fig. [3] By pinpointing the location of the crossing, we 
obtain a first estimation of the ME: E c ~ —0.48. 

Guided by the one-parameter scaling theory of An¬ 
derson localization [22j], we now postulate that Eq. © 
holds not only in the close vicinity of the ME [where 
Xr[E) oc E — E c \ but also away from it, and propose 
to verify this hypothesis by a rigorous finite-size scaling 
analysis of the numerical data for the CBS width. For 
this purpose, we introduce a fitting function of the data 
by Taylor expanding Eq. © up to a certain order tir [i| , 


riR 

A =Y J Xr{E) n L n ' v F n , 


n=0 


( 3 ) 


and further expand the variable Xr{E) according to 
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Xr(E) = X)m=i b m{E c - E) m . Iii this model, F n , b m , 
v and E c are free parameters. We determine them us¬ 
ing a least-square fit of the data for A with Eq. m 
retaining data only for sufficiently long times (such that 
L > 20). We show in the right panel of Fig. [3] the results 
of this fit for curves In A versus E (solid lines). We used 
ur = 2 ,m,R = 3 (that is 7 fitting parameters) for 1141 
data points. The y 2 per degree of freedom is found to be 
0.55. This small value (from the statistical significance 
point of view) comes from the fact that the data got at 
the same energy, but different sizes (i.e. different times), 
are obtained using the same realizations of the disordered 
potential and thus have residual correlations. We have 
also tried to include irrelevant scaling variables to better 
account for deviations to scaling expected at short times 
I,®, but we did not observe significant improvements 
of the quality of the fits. 



FIG. 4. (color online) Left panel: scaling function A con¬ 
structed by fitting the data for A with model ([3]). Points are 
the data, and the solid black curve is the fit. All data lie on 
the same master curve, in agreement with the one-parameter 
scaling hypothesis, Eq. ©■ Right panel: 1 /i(E) = \Xr{E)\ v 
versus energy E (solid curve), together with the confidence 
interval (shadowed region, green online). 1/£(E) vanishes at 
the ME, and is proportional to \E — E c \ v in its vicinity. The 
dashed curve is the prediction obtained from an independent 
finite-size scaling analysis based on the transfer-matrix ap¬ 
proach. 

We then plot the data In A as a function of £(E)/L, 
where i;(E) = \x r {E)\~ l/ (colored points), together with 
the fit to model (j3]) (solid curve). The results are shown 
in the left panel of Fig. Q] We see that all data collapse 
almost perfectly on the same master curve. This result 
demonstrates that the function A, as computed from the 
width of the CBS peak, does follow the one-parameter 
scaling theory, in full agreement with Eq. ©• The quan¬ 
tity £(E) is proportional to the localization length £(1?) 
on the insulating side of the transition, and proportional 
to the inverse of the diffusion coefficient, 1 /D(E), on the 
metallic side. In the right panel of Fig. [4]we show 1 /£(E) 
as a function of energy, as obtained from the fitting pro¬ 
cedure. As expected, 1 /£(E) vanishes at E = E c , which 
signals the divergence of the localization length and the 


vanishing of the diffusion coefficient. The fitting analysis 
also allows us to provide estimations of E c and of the 
critical exponent v. We find E c = —0.4786 ± 13.10 -4 
and v = 1.61 ± 0.03. Because the above-mentioned chi- 
squares are too small, they cannot be used to extract 
the uncertainty. We have thus divided the whole config¬ 
uration sample in several independent subsets, and esti¬ 
mated E c and v for each subset. The reported uncertain¬ 
ties reflect the deviations between the different subsets. 
They are found to weakly depend on tr, most probably 
because the finite size scaling approach relies on data be¬ 
longing to an energy interval much larger than a. In the 
right panel of Fig. U we also display as a dashed curve 
the quantity £,(E) computed from an independent finite- 
size scaling analysis based on the transfer-matrix method 
d, IH]. The latter provides E c = —0.4771 ± 7.10 -4 and 
v = 1.62 ± 0.03, in a somewhat surprisingly good agree¬ 
ment with the estimations extracted from the CBS width. 
The slight discrepancy from the recently reported value 
E c = —0.43 [HI comes from the crude discretization we 
used to save computer resources. Indeed, as involving a 
time propagation and a narrow energy filter, the charac¬ 
terization of the Anderson transition from the CBS peak 
is more numerically demanding than from the transfer- 
matrix approach. Due to this discretization the free- 
space dispersion relation deviates from the massive one 
E = k 2 /2 and p(E) is overestimated near the ME, low¬ 
ering E c . This shift has however no effect on the physics 
of the CBS effect or on the Anderson transition. 

In conclusion, we have shown that the dynamics of 
the CBS peak can be used to characterize the Anderson 
transition, enabling to (i) accurately pinpoint the loca¬ 
tion of the ME (ii) access the critical exponent and (iii) 
test the validity of the single-parameter scaling hypoth¬ 
esis. Our method has the double advantage to be based 
on a physical observable -the CBS peak- which is usually 
well controlled in experiments, and to demonstrate phase 
coherence, which is a crucial requirement prior any claim 
for Anderson localization. The approach has straightfor¬ 
ward applications to the field of atom optics in disordered 
potentials, but it can also be applied to the context of lo¬ 
calization of classical waves [5[ . 
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